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. Abstract. In supersymmetric scenarios with broken i?-parity, baryon number vio- 

(<~^ I lating sparticle decays are possible. We report on the development of a framework 

' allowing detailed studies with special attention given to the hadronization phase. In 

Q^. our model, implemented in the Pythia event generator, the baryon number violat- 

I ing vertex is associated with the appearance of a junction in the colour confinement 

field. This allows us to tell where to look for the extra (anti)baryon directly associ- 
ated with the baryon number violating decay. 
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^ '. 1 Introduction 
O 

CN ' In the Minimal Supersymmetric extension of the Standard Model (MSSM), the standard particle 

^ . content, extended to two Higgs doublets, is doubled up by the presence of superpartners to all 

I normal particles. The conservation of a multiplicative quantum number called i?-parity, defined 

I ' hy R = (_i)2S'+3B+L^ where S is the particle spin, B its baryon number and L its lepton 

number, is usually assumed, since this prevents fast proton decay and has the nice additional 
r~| ! consequence of making the Lightest Supersymmetric Particle (LSP) stable, thus making it a 

WIMP type dark matter candidate. 

However, the choice of i?-parity conservation to prevent fast proton decay is not unique, 
and due to the distinct differences in collider phenomenology between models with and without 
5^ I i?-parity conservation, it is of importance to be well prepared for all possibilities at present as 

well as future high-energy experiments. 

With iZ-parity conserved, experimental SUSY signals would consist of jets, leptons and 
missing E± from escaping neutrinos and LSP's. In scenarios with baryon number violation 
(BNV in the following) the main decay product is jets, with only few leptons or neutrinos, and 
so observability above QCD backgrounds becomes far from trivial at hadron colliders such as 
the Tevatron or the LHC. In order to carry out realistic studies it is therefore necessary to have 
a detailed understanding of the properties of both signal and background events. The prime tool 
for achieving such an understanding is to implement the relevant processes in event generators, 
where simulated events can be studied with all the analysis methods that could be used on the 
real events. 

In this presentation, we concentrate on the possibility that baryon number may be broken, 
resulting in BNV sparticle decays. Sparticle production by BNV, important when the BNV 
couplings are large and/or the sparticles are heavy, is not considered here. In the past, BNV has 
been modelled ^ and studied in detail in the Herwig framework, with emphasis on the 
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perturbative aspects of the production process. In Q, we present a corresponding implemen- 
tation in Pythia, summarized here, where a special effort is dedicated to the non-perturbative 
aspects, allowing us to address the possibility of obtaining a "smoking-gun" evidence that a BNV 
decay has occurred, with questions such as Could the presence of a violated baryon number be 
directly observed? and If so, what strategy should be used?. In addition, many other differences 
exist between the Pythia and Herwig physics scenarios, for parton showers and underlying 
events, thereby allowing useful cross-checks to be carried out and uncertainties to be estimated. 

2 The BNV Scenario 

The most general superpotential which can be written down for the MSSM includes 4 i?-parity 
odd terms: 

•^RPV — 2^zjk(^ ^a^h^ + \jk^ ^a^b + 2 i^l'^ ^ a'i_^a2^a3 + '^t^a^2 Uj 

where i, j, A; run over generations, a,b are SU(2)i isospin indices, and a(j) runs over colours. 

In a S-conserving theory like the SM or the /^-conserving MSSM, there is no colour antisym- 
metric perturbative interaction term, i.e. no term with a colour structure like that of the UDD 
term (the third term in the above equation). Apart from extreme occurrences, like knocking 
two valence quarks out of the same proton in different directions, by two simultaneous but sepa- 
rate interactions, normal high-energy events would therefore not fully display the antisymmetric 
colour structure of the proton. So what is different about the UDD term is that it allows the 
production of three colour carriers at large momentum separation, without the creation of cor- 
responding anticolour carriers. It is the necessary SU(3) gauge connection between these three 
partons that will lead us in the development of the nonperturbative framework. 

A further point about the UDD term is that the contraction of the e tensor with 
implies that A-^-^ should be chosen antisymmetric in its last two indices, since a (j, A;)-symmetric 
part would cancel out. 

The part of the Lagrangian coming from the UDD superpotential term in which we are 
interested is: 

-Cbnv = ^Ar,,6"^"^-3 (^^^^ ^^y^^ (d^)fc + (r (d^)L3 - (j - k)) + h.c. (2) 

where we have made the choice of not yet using any of the antisymmetry requirements, so that 
the ordinary Einstein summation convention applies. 

Combining the vertices in eq. (Q) with the full MSSM Lagrangian, also decays involving one 
or more gauge couplings are clearly possible, e.g. neutralino decay via —>■ qi(— > q^q/cjcii- The 
BNV SUSY decay processes currently implemented in Pythia, with Born level matrix elements 



as calculated by | 


|], are: 




1) djn 


Uidfc 


(36) 


2) iim 


djdk 


(18) 


3) X'n - 


Uidjdk 


(144) 


4) Xt ^ 


UiUjdk 


(30) 


5) Xt - 


didjdk 


(14) 


6) g ^ 


Uidjdk 


(36) 



where n runs over the relevant mass eigenstates: n G {L, R} for the first two generations of 
squarks, n G {1, 2} for the third generation squarks and the charginos, and n G {1, ...,4} for the 
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Figure 1: String drawing in a BNV colour topology. The full lines represent quarks going out 
from the decay vertex, the curly lines gluons emitted in the parton shower, and the dashed 
lines the final strings stretched from each quark across its colour connected gluons back to the 
junction. Note: this picture was drawn in a "pedagogical projection" where distances close to 
the center are greatly exaggerated. 

neutralinos. The numbers in brackets are the number of modes when summed over n, i, j, and 
k, and over charge conjugate modes for the Majorana particles. 

When calculating the partial widths (and hence also the rates) into these channels, we 
integrate these matrix elements over the full phase space with massive b and t quarks, massive 
r leptons, and massive sparticles. All other particles are only treated as massive when checking 
whether the decay is kinematically allowed or not. 

A feature common to the Herwig and Pythia implementations is how double-counting in 
the BNV three-body modes is avoided. The diagrams for these modes contain intermediate 
squarks which may be either on or off the mass shell, depending on the other masses involved 
in the process. If a resonance can be on shell, we risk doing double counting since Pythia is 
then already allowing the process, in the guise of two sequential 1^2 splittings. In particular, 
this means that the list of 1 ^ 3 BNV widths obtained by a call to PYSTAT(2) only represent 
the non-resonant contributions, the resonant ones being accounted for by sequences of 1 — > 2 
splittings in other parts of the code. 

3 BNV Colour Topologies 

Up till now we have considered short-distance processes, where perturbation theory provides 
a valid description in terms of quarks, gluons and other fundamental particles. At longer dis- 
tances, the running of the strong coupling Og leads to confinement and a breakdown of the 
perturbative description of QCD processes. The perhaps most successful and frequently model 
for the transition from the description in terms of quarks and gluons to a description based on 
hadrons is the Lund string fragmentation model Q. 

This approach has not before been applied to the colour topologies encountered in BNV. 
Therefore we here extend the model by the introduction of a junction, where three string pieces 
come together, c.f. figure |l]. Effectively, it is this junction that carries the (anti)baryon number 
that is generated by a BNV process. The hadronization in the region around the junction will 
therefore be of special interest. 

In figure |l|, the central black dot represents such a junction, and the dashed lines show the 
string pieces stretched between the junction and each endpoint quark, across emitted gluons, 
resulting in a Y-shaped topology. In the simplest picture of fragmentation, each string piece is 
broken by the formation of a number of qq pairs along the string. The end-point quark of each 
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piece then pairs up with the closest q (in colour space) to form a meson, leaving a new unpaired 
q which pairs up with another q, and so on until almost all the energy stored in each string 
piece is used up. From this picture, it is evident that the fragmentation eventually produces 3 
unpaired quarks, one on each side of the junction. By colour conservation, with the split off 
mesons being colour singlets, these 3 quarks are in a colour-antisymmetric state, i.e. a baryon. 
In the following, we refer to this baryon as the "junction baryon". 

It could have been interesting to contrast the junction concept with some alternatives, but 
we have been unable to conceive of any realistic such, at least within a stringlike scenario of 
confinement. The closest we come is a V-shape topology, with two string pieces, similar to the 
configuration in a qqg topology. This would be obtained if one e.g. imagined splitting the colour 
(anti-colour) of one of the final state quarks (antiquarks) into two anticolours (colours) . In such 
a scenario the baryon would be produced around this quark, and could be quite high-momentum. 
Of course, such a procedure is arbitrary, since one could equally well pick either of the three 
quarks to be in the privileged position of producing the key baryon. Further, with two string 
pieces now being pulled out from one of the quarks, the net energy stored in the string at a 
given (early) time is larger than in the junction case, meaning the Y junction is energetically 
favoured over the V topology. For these reasons, the V scenario has not been pursued. 

3.1 Fragmentation of Junction Strings 

As mentioned, the kind of string configuration depicted in fig. |^ has not previously been a part 
of Pythia, thus we here outline the technical aspects of the fragmentation process step by step. 
A more comprehensive description will be contained in Q. 

In the rest frame of the junction the opening angle between any pair of quarks is 120°, i.e. 
we have a perfect Mercedes topology. This can be derived from the action of the classical string 
[^], but follows more directly from symmetry arguments. 

Using this requirement, the rest frame of the junction can easily be found for the case of 
three massless quarks (and no further gluons), but the general massive case admits no analytical 
solution. Rather, we use an iterative, numerical procedure. 

When gluon emission is included, the junction motion need not be uniform. Consider e.g. 
an event like the one in fig. |l|. Here the quarks each radiated a gluon, and so the strings to the 
junction are drawn via the respective gluons. It is the direction of these gluons that determines 
the junction motion at early times, and the directions of the quarks themselves are irrelevant. As 
a gluon moves out from the junction origin, it loses energy to the string. From the point when it 
has lost all its energy and onwards, it would then be the direction of the respective quark, and 
not of the gluon, that defines the pull on the junction, resulting in a "jittering around" of the 
junction. Naturally, this also applies in the general case where an arbitrary number of gluons is 
emitted. 

Rather than trying to trace this jitter in detail — which anyway will be at or below the limit 
of what it is quantum mechanically meaningful to speak about — we define an effective pull of 
each string on the junction as if from a single particle with a four-momentum 



Here i = 1 is the innermost gluon, i = 2 is the next-innermost one, and so on up to i = n, the 
endpoint quark. The energy sum in the exponent runs over all gluons inside the one considered 
(meaning it vanishes for i = 1), and is normalized to a free parameter -Enormi which by default 
we associate with the characteristic energy stored in the string at the time of breaking. Note 
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that the energies Ej depend on the choice of frame. A priori, it is the energies in the rest frame 
of the junction which should be used in this sum, yet since these are not known to begin with, 
we employ an iterative procedure. 

Since the string junction is a very localized part of the full string system, it is not desireable 
that the hard part of the fragmentation spectrum of each string, i.e. the hadrons produced 
close to the endpoint quark, should be significantly affected by the presence of the junction. In 
particular, if we consider events where each of the three outgoing quark jets have large energies in 
the junction rest frame, the production of high-momentum particles inside a jet should agree with 
the one of a corresponding jet in an ordinary two-jet event. This can be ensured by performing 
the fragmentation from the outer end of the strings inwards, just like for the ordinary qq string. 
Thus an iterative procedure can be used, whereby the leading q is combined with a newly 
produced qi, to form a meson and leave behind a remainder-jet qi, which is fragmented in its 
turn. Flavour rules, fragmentation functions and handling of gluon-emission-induced kinks on 
the string are identical with the ones of the ordinary string. 

While these hadronization principles as such are clear, and give the bulk of the physics, there 
is a catch: if all three strings are fragmented until only little energy and momentum remain in 
each, and then these remainders are combined to a central baryon, what guarantees that this 
baryon obtains the correct invariant mass it should have? 

In this brief summary, we are forced to refer the reader to Q for the technical details 
pertaining to the answer to this question. The end result is that a physical mass for the junction 
baryon is obtained by first fragmenting two of the three strings from the respective end inwards, 
towards a fictitious other end. In order to have a large-mass system left for the system in which 
energy-momentum conservation will eventually be imposed as a constraint, we prefer to pick 
these two to be the ones with lowest energy, as defined in the junction rest frame. As hadrons 
are successively produced in the fragmentation, their summed energy (in the same frame) is 
updated. Once the hadronic energy exceeds the string energy, the fragmentation has gone too 
far, i.e. it has passed the junction point of the string system, so it is stopped and the latest 
hadron is rejected. 

When two acceptable hadronic chains have been found, the remaining four-momenta from 
the respective two strings are combined into a single parton (diquark), which then replaces the 
junction as endpoint for the third string. If the new parton does not turn out to be spacelike, 
the fragmentation procedure for this string is then identical with that of an ordinary string 
from here on. Otherwise, the fragmentation is restarted from the beginning. Note that popcorn 
baryon production may result in the splitting off of a meson from the initial diquark to produce 
a new diquark. That is, the baryon number may then migrate to higher energies than otherwise, 
but will still be rather centrally produced. 

At this point, it is interesting to see how dependent our model is on the implicit assumptions 
that go into it, for example the definition of the junction pull vector, eq. (^), and the choice 
of the two least energetic string pieces as the ones to be fragmented first in the fragmentation 
scheme described above. 

The variation of the CM momentum spectrum of primary hadrons and junction baryons 
under changes to these assumptions are shown in figure ^ from which it is apparent that the 
model does not suffer from stability problems. Observe also that our earlier remarks that the 
junction baryons would be rather centrally produced are quantified here in the much sharper 
peaking (notice the log scale) of the junction baryon momentum distribution as compared to 
that of the primary hadrons. With respect to the normalization difference between the two sets 
of curves, it is chiefly due to the many mesons produced in the fragmentation. The junction 
baryons in fact roughly double the total number of baryons in the momentum region below ~ 2 
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Figure 2: Momentum spectra of primary hadrons and junction baryons in the decay of a 96 GeV 
neutrahno to three quarks. Results with the default implementation are compared with five 
alternative ones. -Enorm refers to the normalization energy in eq. (|3|). Average multiplicities of 
primary hadrons are shown in the lower right corner of the plot. 



GeV. This gives us our first hint of how to search for this "smoking-gun" evidence of BNV. 

As a final comment, it should be mentioned that more complicated topologies than the ones so 
far mentioned are possible. Specifically when two colour-connected BNV processes occur, there 
will either be two junctions with a string spanned between them or the two baryon numbers 
will cancel against each other and give rise to two unconnected qq string pieces. In the current 
Pythia implementation, we assume that the junction-junction string topology dominates over 
the non-junction one, essentially since we expect the string length, and hence the total string 
energy, to be smaller more often for the former topology than for the latter. 

4 Conclusion 

It has not previously been possible to study baryon number violating decays of SUSY particles 
within the Pythia framework, essentially because it lacked a hadronization mechanism for colour 
configurations containing non-zero baryon number. From Pythia 6.207 on this is now possible, 
and the various aspects of the implementation have been described in broad terms here. Details 
will be available in |^] and in the Pythia manual. The hadronization is based on a physical 
picture and shows negligible model dependence. Furthermore, it allows us to "predict" that the 
smoking-gun evidence of baryon number violation, an excess of baryons, should be looked for in 
baryons having small momenta relative to their parent sparticle. 
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